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Electromagnetic fields of an infinitely long conducting cylinder in a magneto-ionic medium with 
axial static magnetic field are theoretically analyzed. Expressions of the electromagnetic fields and 
the dispersion formulas are obtained. Explicit approximate expressions of the relative propagation 
constants for the case of an extremely thin wire are derived and their dependence on the plasma 
and the cyclotron frequencies is discussed. The behavior of the electromagnetic fields around the 
wire is also discussed. 

1. Introduction 

Extensive investigations have already been made of the electromagnetic wave propagation 
in an anisotropic medium. However, the characteristics of antennas placed in such a medium 
[Kogelnik, 1960; Marini, 1963; Wu, 1963] have not been thoroughly studied. For example, the 
wavelength of the current distribution along an antenna wire has not been obtained, while the 
radiation impedances of linear antennas have already been computed with assumed current 
distributions [Ament et al., 1964]. 

The purpose of this paper is to present some results obtained by a theoretical analysis of 
the electromagnetic waves along an infinitely long and thin conducting wire in a magnetoionic 
medium with the axial static magnetic field [Mushiake, 1964] 

In this theory the electric property of the magneto-ionic medium is assumed to be expressed 
by the well-known tensor dielectric constant, and the existence of ion sheath is neglected. 

Expressions of electromagnetic fields and the dispersion formulas are obtained. Explicit 
approximate expressions of the relative propagation constants for the case of an extremely thin 
wire are derived and their dependence on the plasma frequency and the cyclotron frequency is 
studied. The behavior of the electromagnetic fields of the extremely thin wire will also be 
discussed. 
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Figure 1. Geometry of a circular cylinder. 
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2. Formulation for the Electromagnetic Fields 

The geometry of an infinitely long circular cylinder of perfect conductivity is shown in figure 
1. From Maxwell's equation in a lossless magneto-ionic medium, the following equation for the 
electric field E can be derived: 



where 



VxVx£ = *J(k)-£, (1) 
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X = «>ya>\ (3) 

F=o)///co, (4) 

k 2 = (o 2 iioeo, (5) 

and o>, cup, and o>// are the angular wave frequency, the angular electron plasma frequency, and the 
angular electron gyrofrequency, respectively. The time factor of eft* is assumed, and MKS 
units are used in the present paper. 

The axially symmetric solution of (1) for the electric field can be written as 

E r =EorHf(krr)e-*z z , 

E<t, = Eo<t>Hf\k r r)e-JkzZ, (6) 

E z =E zH<fi\k r r)e-S*zZ, 

where k r and k z are the propagation constants, and the amplitude factors £or, E04,, and Eq z must 
satisfy three relations: 



n 2 z Kz\ 

Ki 4- — E 0r +JK2E06 = 0, 

n 2 r -K 3 / 



(7) 



-JK 2 E 0r + (ki -n 2 r — ntyEot = 0, (8) 

jn z n r E 0r — (n 2 r — Ks)Eq 2 = 0. (9) 

In these equations n r and n z are the relative propagation constants defined by n r = k r /ko and 
n z = k z /ko, and they have complex values of negative imaginary parts in general. Elimination of 
the amplitude factors from (7) and (8) leads to 

Kin*+ (fCi + K 3 )n 2 r n 2 z — (k\-\- K1K3 — Kl)n 2 — 2K 1 K 3 n 2 z J r K3n 4 z -\- (k 2 — k 2 )k 3 = 0. (10) 

This equation gives a relation between the relative propagation constants n r and n z . 

It is found from (10) that two values of n 2 will be determined for one value of n\. In other 
words, two different wave modes which propagate in the same direction of +Z or — Z along the 
cylinder with one common wavelength and which have different propagation constants in r-direc- 
tion can exist. The values of various parameters for these two modes will be discriminated by 
the subscripts 1 and 2 hereafter. By taking two such modes, a general solution of (1) is given by 
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E, = {£on// ( 1 2) (^ir) + £or2// ( 1 2) (A; r2 r)}e^ z , 
E z = {E nH%\k rl r) + E z 2 H%Xkr2r)}e^. 



(11) 



Introduction of the boundary conditions E^ = and E z = into (11) at the surface of the cylin- 
der, or r = a, yields two relations between amplitude factors for two modes. Eliminating the 
amplitude factors from these two relations, we obtain 



ftn(tti— ri\— n^i) H^\k n rl a) _ n^Ki — ri\ — n 2 r2 ) H^XkpTirza) 
' H^(kon rl a) k 3 - n 2 r2 ' H?\k n r2 a) 



K 3 — n 2 rl 



(12) 



(13) 



This is another relation for n r i f ^2, and n z . Therefore, the values of these relative propagation 
constants are definitely prescribed by two simultaneous equations of (10) and (12). This means 
that these equations are the dispersion relation expressed in the form of simultaneous equations. 
Explicit but formal expressions of the electromagnetic fields for this case can easily be obtained 
after some straightforward calculations as follows: 

Er=EAa?Xk«r)- Kl ~ n \~ n ? ■ Snflfiwle^, 
{ Ki — n\— re*, H\ z '(k,-ia) J 

& = E J?— y [mito^TjIsfir* HfXkrAe^*, 

Ki — re| — npj I H\ z \k rZ a) J 

jre r ,re, | //^(^,a) 1 

//,. = - E, ^^ J" W rl r) - ' , , ' WfKkrtr) eW, 

Ki — re|— «£, Vjttot HY>{kr2a) 

H , = £o _S5L. J^ W Arir) - ^§J^ Wrt r)} .*V. 

K l -n z z -n 2 rl V^C n r \H^\k r2 a) J 

where an amplitude factor £0 is to be determined from the total axial current /(z) on the cylinder 
by the relation 



(14) 



/(z) = 27ra[J/*(z)]r=a. 



(15) 



3. Approximations for Extremely Thin Wire 

It seems that the determination of the values of the relative propagation constants from the 
simultaneous equations (10) and (12) is not simple in general. For this reason, only the case of an 
extremely thin wire will be discussed under the assumption of the relations 



\kr\a\ < 1, |A: r2 a| <^ 1, 



1 n rl I ^ I 1 7 I 

In — <^ In k r \a\. 



(16) 
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Under these conditions (12), which is obtained from the boundary conditions, can be approximated 
(see appendix) by 



n 2 x (K\ — n 2 z — h 



n 2 r2 (K l ■ 



ni — n% 



K3 — ni 



K3 — ni 



(17) 



If this approximate equation is used instead of the exact equation (12), the simultaneous equations 
for the relative propagation constants can be solved analytically, and their values are given ex- 
plicitly by 



Kl = l- 



X 



1-Y 2 ' 



n 2 n 



2 r2 



-g-{-K,±VKl + 4KiK 3 } 



XY 



2H-X-Y 2 ) (1-Y 2 ) 



{ -XY± y/X 2 Y i + 4(1 -X) (1 - F 2 ) (1 -X- F 2 )}. 



(18) 



Numerical computations have been made of the relative propagation constants with these 
formulas and some examples of the results are shown in figures 2 and 3. The values of n\ are 
also shown on the plane of X — Y 2 as in figure 4. According to (18), the value of n z is independent 
of the component of the dielectric tensor parallel to the wire, but it obviously depends upon the 
values of X and Y. 

As these results are obtained under the assumption of (16), those assumed conditions must be 
examined. That will be done in the next section. 
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FIGURE 2. Influence of the static magnetic field on the Figure 3. Influence of the static magnetic field on the 

values of the relative propagation constant n z (n\ versus values of the relative propagation constants n r i and 

Y with parameter X). n r2 (nj, and ri z fl versus Y with parameter X). 
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(ii) 
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FIGURE 4. Lines for n 2 = constant on the plane ofX—Y' 1 
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FIGURE 5. Six regions ofX — Y 2 plane. 



4. Discussion of the Relative Propagation Constants 

It is obvious that the assumed relations (16) are not true when either one of \n r i\ > 1, |n r 2| ^ 1, 
\n r i/n r 2\ — or °° holds. For convenience in examining these conditions the X — Y 2 plane is divided 
into six regions as shown in figure 5, where the boundaries of the regions are 



(a) X=l, 

(b) F* = l, 

(c) X+F 2 = l, 

(d) 4(X-1)(F 2 -1)(X + F 2 -1)=X 2 F 2 . J 



(19) 



The approximate values of n 2 . { and n 2 2 in the six regions and on the boundaries are examined from 
the expressions of (18), and the results are listed below, where the value of n 2 is also shown. 



(i) < n 2 x < 1 -X, n 2 2 < 0, < n\ < 1, 

(ii) < n 2 , < 1 -X, n 2 rl < 0, 1< n\, 

(iii) n 2 rl <0, rc 2 2 <0, l<nf, 

(iv) n 2 .j, n 2 2 : complex conjugates, 1 < /if , 

(v) l-X</i 2 1 <l, 0<rc 2 2 , n?<0, 

(vi) ra 2 ,<0, 0<n 2 2 , /i?<0, 
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(20) 
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re|gl. 



(21) 



It is found from these results that the approximate expressions (18) for the relative propa- 
gations constants do not hold on the lines of X=l f Y 2, = \,X J rY 2 = \, and in the vicinity of these 
lines. Further examination reveals that the approximation of (18) does not hold for extremely 
large value of X either. 

Aside from these exceptions, the relations (20) and (21) give various information concerning 
the waves along an extremely thin wire. Some examples are described below. 

When n z is real, the wavelength X along the wire is given, in general, by 



X = XoMz, 



(22) 



where X is the free space wavelength. Therefore, it is found that the wavelength along the wire 
is longer than X in the region (i), while X is shorter in the regions (ii), (iii), and (iv). On the other 
hand, pure imaginary values of n z as in the regions (v) and (vi) represent the attenuating waves in 
the direction of propagation, +Z or — Z. 

If n 2 ri (i = 1 or 2) is positive real, the electromagnetic field components are proportional to 
Hankel functions. On the contrary, if n 2 ri is negative real, the arguments of the Hankel functions 
become pure imaginary and these functions transform into the modified Bessel functions of the 
first kind. Consequently, the field components vanish very rapidly with increasing r. That is 
true, for example, in the region (iii). In the case of the region (i) or (ii), however, the field com- 
ponents corresponding to n r i vanish slowly with increasing r in proportion to the Hankel functions, 
while the components corresponding to n r2 vanish very rapidly with increasing r. 



5. Electromagnetic Fields Near the Wire 



If the conditions 



\k r \r\ < 1, \k r2 r\ < 1 



(23) 



are assumed in the vicinity of an extremely thin wire, electromagnetic fields near to the wire can 
be approximated by 



E r — Eq 



j2(n%-n 2 rl ) e +) k z 



TTnr\(K\ — n 2 z —n 2 n ) k r 



E z — — Zs o 



#r = 0, 



2n r in z | . Th± 
tt(k 3 — nh)\ un. 



[ In Av 2 aJ ' ' 



(24) 
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j2n z K 3 [To n 2 2 - n 2 rl 

7t(k 3 — n^j V/Io ^rl/l,-2 



k r 



H z — Eo 



J2k 2 



7rn r i(Ki — ri\ — n 2 rl ) V ^io 






-h np x In n r \— n 2 r2 \n rc r2 f e+ jA "* z 



(25) 



where 

7 = 1.78107 . . 
Introducing the expression of H& into (15), we obtain 



T/ s r jkn z K 2 (n 2 r2 - n 2 n ) fT . 
Kz) ^ Eo I — { FT V~ e z ' 



(26) 



As the major component of the electric field around the wire is the r-component, the relation 
between £V(r, (/>, Z) and I(Z) might be of interest. For this reason the ratio of E r ioI(Z) is cal- 
culated from (24) and (26). Introduction of (17) into the expression of the ratio finally leads to 



EAr, <f>, Z) _ 1_ 

I(z) 2nrn z 



En. 
e 



(27) 



It should be noticed here that this relation is just the same as in the case of an isotropic medium 
with the relative propagation constant of n z . 

6. Conclusions 

Electromagnetic waves along an infinitely long conducting wire in a plasma medium with 
an axial static magnetic field have been analyzed theoretically. The detailed discussion was con- 
fined to the case of an extremely thin wire. Dependence of the relative propagation constants 
on the plasma and the cyclotron frequencies and their physical significance are explained in 
section 4. The behavior of the electromagnetic fields around the wire is discussed in section 5. 
Although this theory involves approximations and assumptions, some information concerning the 
characteristics of linear antennas in a magneto-ionic medium can be derived from the results of 
this study. For example, the wavelength for the first approximation of the current distribution 
along an antenna wire placed in such a medium may be determined by (22). 



The author acknowledges S. Adachi for his helpful discussions and K. Kato for his assistance 
in the computation. 

7. Appendix 

Derivation of (17). Taking the first terms of the expansions of the Hankel functions in (12) 
under the conditions \k r \a\ < < 1 and \k r 2d\ < < 1 we obtain 



Tl^JKi —Til — n2 n)( K * - ft/V ^ m (Jykr2(ll2) 

■ti 2 2 (ki — ri\ — n 2 . 2 )(K :i — n 2 . x ) In (jyk n al2) 



1 



In (n r iln r2 ) ^_ In (rc r i/ra r2 ) 



In (jykria/2) In k r \a 

Further approximation with the condition |ln (n r i/n r 2)\ < <\ln k,\a\ leads to (17) in the text. 
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